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Abstract. We discuss resonances for Schrodinger operators on metric graphs 
which consists of a finite compact part and a finite number of halflines attached 
to it; the vertex coupling is assumed to be of the 5-type or certain modifications 
of it. Using exterior complex scaling on the graph we show that the resolvent 
and scattering resonances coincide in this case. 



1. Introduction 

Resonances belong to the class of phenomena which are easy to be understood 
on a heuristic level but more difficult when we try to study them rigorously. It is not 
a lack of a precise definition, of course, rather the fact that there are several formal 
ways in which the problem can be approached. The two oldest and most common 
concepts are scattering resonances and resolvent resonances. In the first case one 
inspects the on-shell scattering matrix looking for its sharp changes in some parts 
of the energy axis suggesting a locally enhanced time delay, while in the second one 
we inspect analytical continuation of the resolvent to the "unphysical" sheet (s) and 
look for poles there. The two types of resonances need not be a priori identical, 
and if fact, there is no reason why they should be, because the latter represent a 
property of the Hamiltonian H of the system alone, while the scattering refers to 
a pair (H, Hq) depending thus on the choice of the free operator. Nevertheless, in 
most systems the scattering resonances referring to a "natural" H coincide with 
the resolvent ones, which is a property to be verified in each particular case. 

In this paper we are going to address the question of coincidence between 
scattering and resolvent resonances in noncompact quantum graphs. These systems 
attracted a lot of attention recently from several reasons, the chief ones among them 
being that they can be used to model a wide family of semiconductor and other 
microstructures, and that at the same time they represent a testing ground for 
investigation of fundamental effects such a quantum chaotic behavior; we refer to 
|Ku04,05| for a review and an extensive bibliography to these problems. 
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From the mathematical point of view quantum graphs represent systems of 
ODE's coupled by boundary conditions, so the question posed above can be an- 
swered in a straightforward way, because the resolvent of a graph Hamiltonian can 
be constructed by means of Krein's formula. In reality, however, such a head-on 
approach would be quite laborious because finding the resolvent is not easy unless 
the graphs in question has a trivial topology. Fortunately, there is a simple way: 
using a powerful insight of Jean-Michel Combes and collaborators, first formulated 
in the paper |AC71j . one can transform the search for poles of the analytically 
continued resolvent into the spectral problem for a suitable non-selfadjoint opera- 
tor obtained by an exterior complex scaling of the Hamiltonian. For graphs with 
(5-coupling at the vertices the goal can be be achieved by combing this method with 
the well-known duality property |E97a| . We will first remind some needed notions 
and illustrate our program on examples, then in Section 3 we will formulate and 
prove our main result. 

2. Preliminaries 

2.1. Quantum graphs. The configuration space of our model is a graph with 
a finite number set V = {Xj : j £ 1} of vertices, where / is a corresponding index 
set, and a set of edges attached to them; the graph is supposed to be metric so the 
edges are identified with segments of the real axis. In particular, vertices connected 
to Xj by an edge form the set of neighbors M(Xj) — {X n : n £ v(j) C I\{j}}- 
If necessary we may assume that any two vertices are connected with at most one 
(finite-length) edge, since otherwise we can add a vertex to each "superfluous" 
edge. With this convention, each finite-length is uniquely characterized by a pair 
of vertices, so we can write their set as C = {Cjn ■ {Xj, X n ) £ Ic C I x I}. We will 
assume that a semi- infinite edge is attached to some vertices, denoting the graph and 
its part with edges of finite lengths as T e and T, respectively. The graph boundary 
B is the set of vertices which have a single neighbor, the interior is X = V\B. We 
denote by C the set of vertices to which a semi-infinite link is attached, and we 
introduce Its, It, Ic, respectively, as the appropriate index subsets in /. 

Using the metric structure of our graph we can introduce the state Hilbert 
space of the problem as the I? space on T e which is naturally identified with 

H= L 2 ([O,^])0L 2 ([O,oo)), 

(J,n)elc jeic 

its elements being written by ip = {"4>jn : [j, n ) £ Ic, V'joo : j £ Ic}- The Hamilton- 
ian of the model is a Schrodinger operator on T e which acts as —-T-r+Vj n (x) at each 
edge with a family of potentials {Vj n , Vjoo} and appropriate boundary conditions 
at the vertices. To make things simple we will suppose that Vj n € L°°([0, lj n ]), 
and moreover, that the motion is free at each external link, Vjoo — 0. The describe 
the boundary conditions we need boundary values. Identifying Xj with i = 0wc 
introduce ipj n = lim x _>o+ ipj n (x) and ip'j n = ]im x —>o+ ip'j n (x) , and arrange them as 
into columns, fyj — (ijjj n (Xj) : n £ v(j)) T , and similarly for The bound- 
ary condition at the vertices must ensure that the Schrodinger operator on T e is 
self-adjoint; it is well known IKS!) 9 that this happens if 



(2.1) 
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where Aj, Bj are matrices of dj := card^(j) such that (Aj,Bj) has maximal rank 
AjB* is self-adjoint at each vertex. Moreover, the non-uniqueness in the above 
conditions can be removed if one chooses Aj + Uj — I, Bj = i(Uj + 1) where Uj is a 
dxd unitary matrix HaOO, KSOO . Hence the family of admissible is characterized 
by d 2 real parameters for each vertex. In this large set there are some distinguished 
case. One of them is the 5-coupling for which the wave functions are continuous 
at the vertex, ip 3n {Xj) = ipj m (Xj) =: ipj, and Z)„ £l /(j) ^jni^j) = a j^j- Another 
interesting case is the 5' s -coupling defined by similar conditions with the roles of 
functions and derivatives interchanged; other examples will be mentioned below. 



2.2. Exterior complex scaling. A change of coordinates gives means gener- 
ally to replace the Hamiltonian by a unitarily equivalent operator. The basic idea 
of the method mentioned above is to use a family of such transformation, corre- 
sponding to scalings in the whole space or in the exterior of a given domain, and 
to extend it analytically to complex values of the scaling parameter. This typically 
leads to a non-selfadjoint operator with the essential spectrum rotated; by a suit- 
able parameter choice the resonance poles can be then found as complex eigenvalues 
of the transformed operator, corresponding to eigenfunctions in the Hilbert space. 

For our system a natural idea is to employ a family of transformations which 
leaves the compact part T of T e intact and scales the external semi-infinite edges. 
Identifying such a link with the halfline R + we consider the scaling transformation 
go — > Ugg = e B / 2 g{xe s ) with the parameter 9. Application of Ug to the Laplacian 
leads to its multiplication by e~ 26 , on the halfline one has naturally to transform 
also the boundary condition. Scaling the external edges, all with same parameter, 
we replace therefore the original graph Hamiltonian by the operator 



{fjn}/ \{~fjn + Vjnfjn} 

where the upper component corresponds to the external edges, gj being the wave 
function on the halfline attached to the vertex Xj, and {fj n } corresponds sim- 
ilarly to interior edges of the graph. The domain of the transformed operator 
consists of functions with components fj n S W 22 ([0,lj n ]) and gjg = Uggj with 
gj G W 2 ' 2 (Cj OQ ) satisfying the appropriately transformed boundary conditions. We 
will be interested in the nontrivial situation when 9 is complex, for instance 9 = i"d 
with d > 0. The essential spectrum of Hg comes clearly from the external edges 
and the above formula shows that it is rotates into the lower complex halfplane; for 
i9 large enough one can uncover the poles of H laying on the second sheet. 

Before we pass to formulating the general result, let us discuss several examples. 



2.3. Example: a line with an appendix. The simplest nontrivial example 
consist a line, representing to external links, to which a line segment of length I > 
is attached at the point x = 0. Consequently, the Hilbert space is L 2 (R) ©L 2 ([0, 1]), 
its elements can be written as -0 = (?) where g, f refers to the line and the appendix, 

respectively, and the Hamiltonian acts as Hip = (_ r7/?j_y t) ■ To make it self-adjoint 
we have to choose boundary conditions at the point x = and the other endpoint 
of the segment. Let the latter be Dirichlet, f(l) = 0, and at the vertex we put 



g(0+) = g(0~) =: 3(0), /(0) - 0g(O) + 7 /'(0), </(0+)-</(Q-) = Sg(0) — Af'(G) 



4 



PAVEL EXNER AND JIM LIPOVSKY 



for some /3,7,<5 6 K following ES94) : it is the most general class of time-reversal 



Hamiltonians with the line wave function continuous at x = 0, in particular, the 
case /? = 1, 7 = represents t he ^-co upling with a — 5. The scattering problem 
for this system is easy to solve |ES94| : denoting by fi(x) the normalized Dirichlet 
solution to the Schrodinger equation at the appendix, we can express the reflection 
and transmission amplitudes at energy k 2 as 

<5[/ ; (0)- 7 //(0)]-/3 2 //(0) t= 2t*[/ I (0)- 7 //(0)] 



(2*fc - 8)[fi(0) - 7 //(0)] + /32//(0) ' {2ik - 5)[fi(0) - 7 //(0)] + ^//(0) ' 

Scattering resonances coincides with zero s of the denominator; their behavior with 
respect to the parameters is discussed in ES94] . If we suppose for simplicity that 
V = 0, then fi(x) = sin k(l — x) and the condition can be rewritten as 

3 2 k 

(2.2) tanW=— — ~fk. 

2ik — 5 

Let us look now at the problem from the complex scaling point of view and put 
ge{x) := e 9 / 2 g{e 9 x). The corresponding boundary values are ge(0) = e e ^ 2 g(0) 
and g' e {0±) — e 36 / 2 g'(0±). This has to be substituted into the above boundary 
conditions: the continuity at x = and the Dirichlet condition at the appendix 
endpoint do not change, while the other two yield 

/(0) = (3e- 9 / 2 g g (0) + 7 /'(0) , e- 39 / 2 [g' 9 (0+) - g' e (0-)} = Sc' 9 / 2 (0)g e (0) - /3/'(0) ; 

this can be regarded as the boundary conditions which define the non-sclfadjoint 
operator Hg. Assuming again V = 0, it is easy to solve the corresponding spectral 
problem. The appendix solution is f(x) — b sinfc(/ — x) and its halfline counter- 
parts are the corresponding exponential functions e^ lkxc . The first one of the 
above conditions gives 

Be- 9 / 2 

and substituting from here into the second one we get 

l3 2 (-kcoskl) 



- 39 / 2 ge(0)2tke 9 =g e (0)e- 9 / 2 



sin kl + 7fccos kl 



It is straightforward to conclude from the last relation that the resolvent resonances 
for the line with an appendix are determined by the condition (|2.2(l again. 

2.4. Example: a loop with two leads. As the next example let us consider 
a graph consisting of two internal edges of lengths Zi, Z 2 , connecting the endpoints 
of two halflines. Consequently, the Hilbert space is L 2 {Mr) ©L 2 (R+) L 2 ([0, h])® 
L 2 {[0, h]) and states are described by the columns ip — (/, g, u, v) T , with the Hamil- 
tonian acting as ip = (— /", —g", —u", —v") T ; for simplicity wc suppose again that 
the particle is free away of the vertices. The boundary conditions will be now chosen 
as a (5-coupling at each vertex, i.e. the continuity together with 

u'(0) + v'(Q) - /'(0) = af(Q) , -u'{h) - v'(l 2 ) + g'(0) = pg(0) . 

The scattering problem is again easy to solve. One uses the above boundary con- 
ditions to match the solutions e lkx + re _ltl and te lkx for ^fx > with the linear 
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combinations of the exponentials at the internal links; solving the corresponding 
system of linear equations we get the on-shell reflection and transmission amplitudes 



i - 7 (fc) 2i 
r ( k ) = —, — 77T' H fc ) " 



where 

7 (fe) := 



i + 7(fe) ' 7(fc) + i 

( 1 , l ^ 2 

I sin/c^i sinfc/2 / / 1 1 \ Ck 



^— i Vtan/cZi tanfc/2/ 



tan/c/i tan fc^2 

In particular, the scattering resonances are now determined by the condition 
(2.3) 7(fc)+z = 0. 

Let us now apply the scaling to both halflines putting fg(x) := e e ^ 2 f(xe e ) and 
geix) := e e / 2 g(xe e ). Calculating the new boundary values we can write down the 
boundary conditions determining Hg, the continuity at the vertices together with 

-e- 39 / 2 /«(0-) + u'(0) + v'(0) = ae- e / 2 f e (0-), 

e- 36/2 g' e (0+) - u%) - v'(l 2 ) = (3e- e / 2 g 9 (0+). 

To solve the eigenvalue problem for Hg one has to match the solutions at the 
external edges which are e Tlkxc with 

_ /g(O-) sinfc(/i - x) + gg(0+) sinfcx _ e/2 
u(x) — : — e 

sin fcti 

and the analogous solution at the other internal link. This yields the conditions 

= (a-ik)f e (0-) , 



gg(0+) [ -\t- + -^—) - / fl (0-) \-\r + -\r 
1 sinfeti smfc/2/ \tanfc/i tanre<2 



(p-ik)g e (0+) , 



t&nkli tanfc^/ sinfc^ 
which are obviously equivalent to the scattering- resonance equation l|2.3[l . 

2.5. Example: a magnetic lasso graph. To show that these consideration 
extend beyond the pure Schrodinger case, take now a graph consisting of a loop 
of circumference / to which a halfline is attached and suppose that it is placed 
into a magnetic field perpendicular to the loop plane. The Hilbert space is thus 
L 2 (R + ) (J) L 2 ([0, 1}) with elements ip = (j) and the Hamiltonian acts as 

H * = H { 9 i) = (-/»- Z r + a 2 L 

where A is the corresponding vector potential. In fact, the form of the field is not 
important, what matters is the flux through the loop. The boundary conditions we 
consider are similar to those of Example 12. 31 wc suppose that the wave function is 
continuous on the loop, f(0) — /(/), and furthermore, 

/(0) = crVW - /'(/)] + 7S'(0) , g(0) = 7[/'(0) - f'(l)} + cTY(O) 

for some a, a £ R and 7 £ C; now we do not require time-reversal invariance. 
Performing the complex scaling on the halfline in the same way as above, the 
boundary conditions become 

/(Q) = a-Vfl)) -f(l)] +1 e- 3 ^ 2 g' e (0), 

e- 9/2 ge(0) = 7[/'(0) -/'(/)] +5- 1 c- 3 ^(0). 
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To solve the eigenvalue problem for Hg we use the following Ansatz: ge{x) — e zkc 
and f(x) = Ce~ lAx sin (kx + ip) with tp given by tamp — sinkl(e lAx — coskl). A 
straightforward computation yields then the resonance condition 

fk ik 2 \i\ 2 \ 

(2.4) sinfcZ-2 - + - ^— r (cos - cosfcZ) =0, 

\ct 1 — ika 1 J 

which is again the same as the following from the scattering on the lasso |E97b| . 

3. General graphs 

3.1. Choice of the vertex coupling. In order to extend the observation 
made in the examples to the class of graphs described in Sec. 12.11 we have first to 
specify the vertex boundary conditions we will consider, in particular, the way in 
which the external links are attached to T. For simplicity we suppose that at most 
one halfline sprouts of each point of the boundary of T. Let us denote the wave 
function on such a halfline referring to Xj € B as gj, and those on edges joining Xj 
with interior vertices as fj n ,j = l,...,m. To keep things simple we will consider 
a coupling which generalizes directly the one of Example 12.51 the functions are 
continuous at the vertex, fji(0) = /ia(0) = • • • = /jm(0) ='■ /j(0), and 

m m 
(3.1) fj(0) = aj 1 f jn (°) + 7jflK°) ' 5,(0) = 7, ]T f' jn (0) + afg^Q) 

n—l n—1 

for aj, ctj € R and jj € C; it is straightforward to check that it is of the type 1)2.1(1 . 
In the interior vertices of T we suppose a S coupling; in the same notation it is 
continuity again, f 3 i(0) = /j2(0) = . . . = fj m (0) —■ fj(0), and 

m 

(3-2) E^(°) = a i/i(°) 

n=l 

for a real aj, in general different at different vertices. Finally, if there is a vertex 
in the boundary of T to which no external edge is attached, we assume Dirichlct 
boundary conditions there, fj n — for j G Ib\Ic- 

3.2. A duality. As we have indicated in the introduction, our second main 
tool will be a duality between Schrodinger equation, 

(3.3) Hip = k 2 ip , 

on a graph with a (^-coupling and a certain difference equation. Let us recall it 
now in more details. To cover both T and T e it is useful to formulate the result 
in a way which allows to describe generalized eigenfunctions at the same time, 
hence we consider the class D\ oc (H) which is the subset in \J ^ n \ e j c L 2 (0, ijn) (the 
direct sum) consisting of the functions which satisfy all the requirements imposed 
at ip £ D(H) except the global square integrability. 

On the edge C n j = [0,£j n ], with the right endpoint identified to the vertex Xj, 
we denote as Uj n and Vj n the normalized Dirichlet solutions to — /" + Uj n f — k 2 f, 
i.e. those satisfying the boundary conditions 

U jn (£jn) = l-{ U j n )'{lj n ) = , V jn (0) = l-(vj n )'(0) = ; 

their Wronskian is Wj n = —Vj n (£jn) = Uj n (0). Since in our case the set V is finite, 
the assumptions used in |E97a| are satisfied and we have the following result. 
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Theorem 3.1. Suppose that ijj G Di oc (H) solves for some k ft IC with 

k 2 G K, Imk > 0. Then the corresponding boundary values satisfy the equation 

(34) v ^.(^ia^l^^.o, 



and conversely, any solution {tpj : j G Ix} of the system to \'d-4\l determines a 
solution of by 

Tpjnix) = tt7~~ u jn(x) ~ V jn (x) if n G u(j) D Jz , 

ib ■ 

tpjn(x) = - —y- vj n (x) if n S n I e ■ 

When applied to r e the above result concerns also generalized eigenfunctions. 
It is useful, however, to specify it for the scattering situation. Let us consider the 
solutions on the external links gj(x) — a,je~ lkx + bje lkx for all j 6 Ic- The operator 
S maps the vector of incoming amplitudes a = {<Xj} into the vector of outgoing 
amplitudes b = {bj}, i.e. b = Sa. Poles of the scattering matrix are given by the 
condition det S 1-1 = 0; recall that for a nonreal k the matrix ceases to be unitary. 
Substituting the Ansatz into the boundary conditions H3.ll) we obtain for j G Ic 



n=l 

m 

& j( a j + b j) = & j7j Y ^jn(°) + ik ( b j - aj) , 



n=l 

while for j fi Ic we have the standard 5-coupling 

a j4>jU) = Y i>'jnU) ■ 
nev(j) 

Using Theorem [O] we can now proceed in the way similar to |E97b| to obtain the 
system of equations for j 6 Ic 

(3.5) a s 1>j= Y + H ^ ^ + ''•'"■<' ' Ul J ~ a o) » 

/ T *' in . .. yy jn 



(3.6) ajibj + aj ) = a j7j ( Y ~W~ + Y Vj W^ ^ I + ik ( b i ~ a ^ ■ 

On the other hand, for j ft Ic we have the condition (|3.4() . We have thus arrived 
at a system of N = / + Ic equations for ibj and bj which gives, in particular, the 
sought S-matrix relating the incoming and outgoing amplitudes. 

3.3. Duality for a complex-scaled graph. Let us now perform the exte- 
rior complex scaling on T e which changes the external-edge wave function gj to 
9j,e{ x ) '■= ed ' 2 9j( e0 %)i the scaling parameter will be the same for all j G Ic- Our 



s 
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aim is to find eigenvalues of the scaled operator He, and since there are no poten- 
tials on the external links, we know that the solutions there will be of the form 
9j,e(x) = e lkc x gj,e{0)- Substituting from here to (|3.1(1 we get 

m m 

/i(0) = oq 1 £ /j»(0) + ^e- e / 2 9j , e (0) , e^g,, e (0) - 7j £ /j„(0) + iter 1 , 

n— 1 n— 1 

and eliminating #7,0(0) from here we arrive at 

«»>=(v +T ^)xx.(°>- 

\ J / n=l 

In other words, the scaling to replacement of the coupling at the vertex by a new 
efficient one, non-self adjoint and energy-dependent, with the parameter 

1 — zkot ■ ^ 

^ 37 ' ) ^ k aJ l + ifc(| 7i | 2 Q J -6q l ) ' 

which applies only to the interior edges meeting at the vertex Xj. Now we can 
repeat step by step the proof of Theorem 13. II given in E97a] to find the system of 
difference equations determining the eigenfunctions of Hg through their values at 
the vertices. Those referring to the interior ones, j 6 Ij, do not change being again 
given by (|3.4|) . On the other hand, the equations for j 6 Iq become 

(^) E r-f^ %^-/^)U = o, 

ne»U)nix Wjn \neu(j) W > n ) 

Now we are ready to compare both systems. Substituting from 1)3. 5JI into l|3. 6J1 we 
get after a straightforward computation 

b _ a <^rwh - 2oj 

3 3 l + iHa^-aJ 1 )- 
Substituting this into 13.5(1 again we obtain the system of equations 

y y ^-ftJi^ 2 lkam 

$t T W jn \ ^ W jn MJIfj 1 + ife (| 7 |2 

and 1(3. 4JI . As a final step, it is easy to see that the determinants of both system 
yield the same pole condition. 

3.4. The main result. Summarizing the above discussion we are able now 
to state the claim announced in the introduction about the two resonance sets. 



Theorem 3.2. Let H be a Schrddinger operator on T e as described in Secs. \2.l\ 
and \3.1\ then the families of its resolvent and scattering resonances coincide. 
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4. Concluding remarks 

The situation when the graph loops are pierced by magnetic fluxes the example 
of which we saw in Sec. 12. 51 can treated in the general case also since it reduces to a 
simple transformation of wave functions on the internal edges E97a . Furthermore, 
the duality used in the argument is not restricted to the (5-coupling; in E97a its 
validity is demonstrated for its ^-counterpart. On the other hand, it is natural 
to expect that the equivalence of the two resonance types is valid for any coupling 
(|2.1|) . however, to prove this claim the present approach needs to be modified. 
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